Sliding phase in randomly stacked 2D superfluids/superconductors 
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Using large scale quantum Monte Carlo simulations of lattice bosonic models, we precisely 
investigate the effect of weak Josephson tunneling between 2D superfluid or superconducting 
layers. In the clean case, the Kosterlitz-Thouless transition immediately turns into 3DXY, 
with phase coherence and superflow in all spatial directions, and a strong enhancement of the 
critical temperature. However, when disorder is present, rare regions fluctuations can lead to an 
intermediate finite temperature phase — the so called sliding regime — where only 2D superflow 
occurs within the layers without any transverse superfluid coherence, while a true 3D Bose-Einstein 
condensate exists. Critical properties of such an unconventional regime are carefully investigated. 



Introduction — In condensed matter quantum systems, 
disorder is known to lead to a large variety of fascinating 
phenomena. Relevant in strictly one-dimensional (ID) 
systems [1], the situation is usually less universal in 
higher dimension [2]. In that respect, the case of 
disordered superfluids or superconductors has been 
intensively studied during the past decades [3]. More 
recently, ultracold atom experiments have triggered 
increasing interest for disordered induced localized 
phases of interacting bosonic systems [4]. In such a 
context, the question of the competition between super- 
fluidity or superconductivity and disorder has led to a 
large amount of fascinating works focusing on bosonic 
quantum glasses. In particular, Bose glass physics has 
been addressed in various condensed matter systems: 
superconductors, cold atoms, polaritons, quantum anti- 
ferromagnets [5]. Among all these recent developments, 
it appears that the subtle interplay between disorder 
and dimensional crossover has been overlooked, despite 
the strong experimental interest in layered systems 
such as superconducting cuprates [6], thin superfluid 
films [7], or quasi-2D bosonic gases [8]. Only recently, 
two simultaneous analytical works [9, 10] have predicted 
for weakly disordered coupled superfluid layers an 
intermediate sliding regime (first evoked in the context 
of DNA complexes [11]), where there is no phase locking 
between 2D superfluids. Such exotic sliding superfluid 
state had also been pursued in the context of the frus- 
trated antiferromagnetic layered system BaCuSi206 [12]. 
However, frustration-induced 2D classical decoupling 
does not survive quantum fluctuations which restore 
3D superfluidity [13, 14]. Dynamical decoupling of 
superconducting layers induced by a stripe order has 
also been discussed for I^-^Ba^CuCU [15]. 

In this work we focus on weakly disordered superfluid 
or superconducting layers. To this aim, we use large 
scale quantum Monte Carlo (QMC) simulations of lat- 
tice bosonic models and test predictions of Refs. [9, 10] 
regarding the occurrence of a stable sliding regime for 
disordered coupled superfluid layers. We first present 



evidences that when isolated layers are coupled by a 
Josephson tunneling t±, a finite transverse stiffness im- 
mediately develops for any t± ^ 0. On the other hand, 
when layers of two types, having different individual 
Kosterlitz-Thouless temperatures, are randomly stacked 
and 3D coupled, QMC simulations show evidences for 
a wide temperature window where a sliding phase is 
achieved. We then present a detailed study of the critical 
properties of such a state, and discuss momentum space 
properties. 

Coupled clean superfluid layers — One of the simplest 
lattice model for studying 2D superfluidity is the Bose- 
Hubbard Hamiltonian on a L x L square lattice 

- /j J ^2n i ^U^2n i (n i -l). (1) 

i i 

Particle hopping is restricted to pairs (ij) of neighboring 
sites, and the density of bosons (n) = Yli^ib^/L 2 
is controlled by the chemical potential fi. In the 
following, we will work in the limit of infinite on-site 
repulsion U — » oo: the hard-core bosons limit, know as 
the Tonk- Girardeau regime, as achieved in cold atom 
experiments [16]. Ground-state and finite temperature 
properties of 2D hard-core bosons on a square lattice 
are pretty well-known [17]. For \fi\ < 2ty, the system 
is superfluid and Bose-condensed at T = whereas 
at finite temperature, only superfluidity survives [18] 
below a finite Kosterlitz-Thouless (KT) temperature 
Tkt where the superfluid density p s f displays a universal 
jump [19]. 

When superfluid layers are Josephson coupled by 

^i = -yE ( b tAe+i + b[e + Ae) > (2) 

where t stands for the layer index, a true 3D long-ranged- 
ordered phase is expected at finite temperature, where 
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FIG. 1: (color online) QMC results for the 3DXY transition of clean coupled superfluid or superconducting layers of hard-core 
bosons [Eqs. 1 and 2]. In-plane and out-of-plane superfluid stiffnesses are shown in panels (a-b) for three representative values 
of t±/t\\, and different system sizes: L = 8 (circle), L = 16 (square), L = 24 (diamond), L = 32 (up triangle), L = 
48 left triangle). In panels (c-d), the same QMC data multiplied by L display a crossing at the critical temperature T c . 



the KT superfluid-only regime immediately turns into 
a Bose condensed superfluid with 3D phase coherence. 
The 3D critical temperature can be estimated from rather 
simple arguments: If one assumes decoupled superfluid 
layers, approaching the KT transition from above the 
correlation length rapidly increases as 



C 2D (T)^Coexp^^-j. (3) 

A crossover to the 3D regime is expected when £_l£|d — 
ty, i.e. below 

r »=^( i+ ,^iw)- (4) 

Interestingly, such an estimate turns out to give the cor- 
rect behavior for the actual T c at small tj_ , as visible 
in Fig. 2 where QMC results are displayed for weakly 



coupled 2D superfluid layers of hard-core bosons at half- 
filling (/i = 0). These results, shown in details in Fig. 1, 
have been obtained using the SSE algorithm [20] at finite 
temperature for various 3D arrays of N = L 3 sites with 
L = 8, • • • ,48. Taking an intralayer hopping strength 
tn = 1, several orders of magnitudes for inter layer Joseph- 
son couplings have been explored: t± — 10 -3 , • • • ,1. For 
weak amplitudes t±/t\\ < 1/10, the critical temperature 
obeys the simple prediction Eq. 4 with Tkt = 0.342t||, 
b = 1.86 and = 0.34, as displayed in Fig. 2. It is worth 
noting the strong logarithmic enhancement of the tran- 
sition temperature above Tkt (for numerical values, see 
table I). 

The superfluid response along all spatial directions 
can be probed by imposing twisted boundary condi- 
tions along a given direction, which translates into the 
computation of winding number fluctuations in the 
QMC scheme [21]. The in-plane superfluid stiffness pi 
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TABLE I: Numerical estimates from Fig. 2 for the 3DXY transition temperature T c /t\\ of clean superfluid layers as a function 
of the interlayer hopping t±/t\\. 



is plotted together with the transverse response pjr in 
the panels (a) and (b) of Fig. 1 for t±/t\\ = 1, 0.1, 0.01. 
There, one clearly sees that quantum coherence of the 
superfluid establishes in both directions, even for very 
weak t±/t\\. When normalized with respect to tn or 
the superfluid response displays a small anisotropy, only 
observed for infinitesimal t±. Importantly, longitudinal 
and transverse stiffnesses vanish at the same critical 
temperature T c . This is best visible in Fig. 1 (c-d), where 
the 3DXY transition is detected by using the critical 
scaling of the stiffness p s ~ L 2 ~ D ~ Z where D = 3 and 
z = 0. Indeed, QMC data of p s x L display nice crossing 
features for different system sizes for both components: 
longitudinal pi and transverse p^. 




FIG. 2: (color online) Behavior of the 3DXY transition tem- 
perature T c /t\\ of clean superluid layers, plotted against inter- 
layer hopping t±/t\\. Green crosses are the QMC estimates 
obtained from crossing analysis displayed in Fig. 1 and the 
red line is a fit to the form Eq. 4 using T K t = 0.342ty [19] 
(green star), b = 1.86, and £o = 0.34. 



Randomly stacked layers: sliding regime — As pre- 
dicted in Refs. [9, 10], an unconventional intermediate 
temperature regime with 2D-only superfluidity can be 
expected in a disordered layered system, if disorder is 
introduced in a way that rare regions have distinct KT 
temperatures. In order to investigate precisely such a 
physics, we focus on a simple 3D quantum model of hard- 
core bosons where L layers of two types (A and B) are 
randomly stacked and weakly coupled: 

L ^ 

e=i (ij) i 
+ |E( 6 lA,i + i + 4 + i^)]> ( 5 ) 

i 

where the in-plane hopping tjj = 1 (A) or 2 (B) with 
probability 1/2, and a constant interlayer Josephson tun- 
neling t± = ty/lO. In the following, we fix the chemical 
potential p = such that the system remains at half- 
filling. Taken independently, each layers exhibit, when 
decoupled, individual Kosterlitz-Thouless temperatures 
T^ T = Tj5? T /2. When the transverse coupling t± is 
turned on, rare thick slabs of M consecutive layers of 
the same type (A or B) appear with a probability ' . 
In the infinite system size limit, the existence of such rare 
regions define an upper temperature 

T u = T| T [1 + 4&7ln 2 (tL$/t)] ~ 1.06i||, (6) 

and a lower temperature T\ =T u /2. In this temperature 
range, QMC simulations (Fig. 3) show clear evidences 
for an intermediate phase T s < T < T u where the 
longitudinal superfluid response p\ is finite whereas the 
transverse one pjr = 0. Numerical results have been 
obtained on very large cubic clusters of size L x L x L, 
with L = 16, 24, 32, 48, 64, 100, and averaged over sev- 
eral hundreds of independent disordered configurations 
of random stacking. Contrary to the clean situation 
(Fig. 1) where, even for tiny transverse tunneling 
a full 3D coherence was found below a single ordering 
temperature T c , with a finite superfluid stiffness in 
all spatial directions, in the random stacking case the 
transverse superfluidity becomes non-zero only for 
T < T s ~ 0.7tn. Below T s the finite value of pjr signals 
a phase locking between layers. In the present situation, 
a lower bound for T s (better than T/) is T^ T ^ 
as given by rare events where the effective tunneling 
between neighboring B-layers — » 0. On the other hand, 
close to T n , infinite size extrapolations of pl(T) are very 
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FIG. 3: (color online) QMC results for the superfluid responses in the longitudinal pi (Left) and the transverse p^ (Right) 
directions plotted against the temperature T/t\\ for different 3D systems L x L x L of size L — 16, 24, 32, 48, 64, 100. Disorder 
averages have been performed from more than 10 4 samples for L — 16 to 10 2 samples for the largest lattices, such that error 
bars are typically smaller than symbols sizes. When extrapolated to the infinite size limit (yellow symbols), the superfluid 
stiffness in the plane direction pi vanishes for T c = 1.04(2)£|| whereas the response in the transverse direction pjr has already 
vanished for T s ~ 0.7£||, thus defining a sliding regime (shaded area) for 0.7 <T< 1.04(2). In the right panel, the transverse 
stiffness for disorder-free A layers with t±/t\\ =1/10 are also shown, defining Ti (see text). 



well described by the analytical prediction of Ref. [9] 
- exp [-a(T c - T)~ v ) with a = 0.6(1), T c ~ 1.04(2) and 
i/ = 0.6(1). 

Scaling of the transverse stiffness — In the sliding 
regime T s < T < T c , the absence of transverse superfluid 
response leads to the following finite size scaling [9, 22] 



H _ 1 &>E{ V± ) 



cp ± =0 



ocL 1 -^). (7) 



Here cp± is a twist angle imposed in the perpendicular 
direction, and z±(T) is the anomalous exponent for 
transverse superfluidity. Finite size scaling for p^-(L,T) 
is displayed in the inset of Fig. 4 for various temperatures 
in the range 0.60 < T/tn < 0.98. In the main panel of 
Fig. 4, estimates for the anomalous dynamical exponent 
z±(T) are plotted and display a clear divergence when 
T approaches T u . 

A precise fit is not reasonably possible since both 
power-law and exponential divergences describe correctly 
the data on such a reduced scale. From Refs. [9, 22], we 
expect a power-law for SU(2) symmetry and exponential 
in the present U(l) case. 

Momentum space properties — An interesting issue lies 
in the momentum space properties. A first observation 
concerns the relationship between Bose-Einstein conden- 
sation and superfluidity. Although the relative phases 
can fluctuate from one layer to the other in the sliding 
regime, they remain locked in each plane. Therefore, a 



macroscopic occupation of the k = mode is expected, 
as visible in Fig. 5 where the condensate fraction 



(8) 



i,3 



is displayed together with the superfluid fraction 
(directly obtained from the average over all spatial 
directions of the superfluid stiffnesses [23]). From a 
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FIG. 4: (color online) Anomalous exponent z±(T) of the finite 
size scaling of the transverse stiffness Eq. 7 estimated from 
QMC data shown in the inset. The blue line is a fit to a 
power-law oc \T — T u \~ a with a = 0.94 and the dashed red 
line is an exponential fit oc exp(a|T — T u \~ u ) with a ~ 0.33 
and v ~ 0.67. 
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FIG. 5: (color online) QMC results for SF and BEC fractions 
of the 3D disordered layered system. Infinite size extrapola- 
tions are shown (yellow symbols) for both fractions. Inset: 
ratio between L — » oo SF and BEC fractions. The horizontal 
line is p s f/po = 1 and the dark curve is a guide to the eyes. 



macroscopic point of view, below T u the disordered 
layered system displays both superfluidity and Bose- 
condensation, without clear evidence for an absence 
of transverse superfluid response. Nevertheless, inside 
the sliding regime there is an anomalously small BEC 
fraction as compared to the SF fraction (see inset of 
Fig. 5), whereas below T s where full 3D coherence 
is recovered, both fractions are of the same order of 
magnitude, as expected for a conventional superfluid. 

The investigation of the momentum space distribution 
JV(k) = 5>t 6r ,)exp[ik.(r-r')] (9) 

turns out to be even more instructive. In Fig. 6, QMC 
results for the distribution N(k) of a single disordered 
sample (a random stacking of 100 2D layers of size 
10 x 10) are plotted against k±, the wave vector in the 
transverse direction, for 3 representative temperature 
regimes: (a) normal phase, (b) sliding regime, (c) 3D 
superfluid. Whereas there is no qualitative difference 
between clean and disordered systems for (a) and (c), 
the sliding regime (which is similar to (c) in the clean 
case) features a much broader momentum distribution in 
the disordered case. We believe that such a qualitative 
feature would be observable using time of flight imaging 
in cold atom experiments. 



Conclusions — Before concluding, we may briefly discuss 
the generality of our results. A schematic tentative 
phase diagram is given in Fig. 7 in the t±/t^ — T/t^ 
plane for model (5). The key ingredient to observe the 
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FIG. 6: (color online) QMC results for the momentum distri- 
bution functions of a single random stacking of 100 A and B 
layers (see inset a) of 10 x 10 sites, plotted against the trans- 
verse momenta k± for various longitudinal wave vectors k\\, 
as shown on the plot. 3 representative temperature regimes 
are displayed for the disordered system (o) and are compared 
with results for a clean sample (•) with tu —2. 



sliding regime, valid for other kinds of random stacking 
or disorder, lies in the existence of a finite temperature 
window T s < T < T u where exponentially rare regions of 
infinite extend having the smallest ordering temperature 
remain disordered whereas other layers are superfluid. 

To conclude, we have numerically studied the effect 
of a weak 3D Josephson tunneling between superfluid 
or superconducting layers. Contrary to the case of 
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FIG. 7: (color online) Schematic phase diagram for randomly 
stacked superfluid layers of type A or B [Eq. (5)]. 

clean identical planes, when distinct layers are randomly 
stacked an intermediate sliding regime emerges over a 
finite temperature window, as predicted recently [9, 10]. 
The absence of 3D phase locking is accompanied with 
a diverging exponent z±(T —> T c ), a finite but very 
small BEC fraction, and a broadening of the momentum 
distribution iV(k). 

I thank F. Mila and T. Vojta for usefull comments and 
discussions. 
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